Gap Probabilities for Double Intervals in Hermitian Random Matrix 
Ensembles as r-Functions Spectrum Singularity case 
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The probability for the exclusion of eigenvalues from an interval 
(—a;, x) symmetrical about the origin for a scaled ensemble of Her- 
mitian random matrices, where the Fredholm kernel is a type of 
Bessel kernel with parameter a (a generalisation of the sine ker- 
nel in the bulk scaling case), is considered. It is shown that this 
f^) ■ probability is the square of a r-function, in the sense of Okamoto, 

C") | for the Painleve system Pm . This then leads to a factorisation of 

the probability as the product of two r-functions for the Painleve 
system Pup . A previous study has given a formula of this type but 
^f-) ' involving Pup systems with different parameters consequently im- 

\Q \ plying an identity between products of r-functions or equivalcntly 

■ sums of Hamiltonians. 

. The probability £73(0; J; g(x); N) that a subset of the real line J is free of eigen- 

(f-) ' values for an ensemble of N x N random matrices with eigenvalue probability 

density function proportional to 

AT 



> 



J: 



Oh- (i) n^o n k--^ 



x 



i=l \<j<k<N 

(j3 = 1, 2 or 4 according to the ensemble exhibiting orthogonal, unitary or symplectic 
symmetry respectively) is a fundamental statistic in the study of these ensembles. 
Most effort has focused on the case where J is a single interval, one endpoint fixed at 
the edge of the support of the measure defining the ensemble whilst the other is free, 
and taken to be the independent variable in the system of equations determining 
the gap probability. There is, however, another interesting case where the set J 
consists of two disconnected intervals, but related to each other so that there is 
still only one free variable. For instance there is the result for unitary ensembles of 
Hermitian matrices that the gap probability for an interval symmetrical about the 
origin J and an even weight function g 2 (x) (where the integrable examples include 
Gaussian, symmetric Jacobi and Cauchy weights) factorises ^ 

(2) E 2 (0;J;g 2 (x);N) 

= E 2 (0; J + ;y-^ 2 g2(y 1/2 ); l(N+l)/2\ )E 2 (0; J+; y 1 ' 2 g^y 1 ' 2 ), \N/2\) , 

( J + is the positive member of a pair of intervals composing J). Examples of where 
this relation has been useful can be found in the above reference and |S] . An example 
of such a double interval statistic is one arising from an ensemble of random unitary 
N x N matrices with the joint eigenvalue pdf 

N 

(3) p(z 1 ,...,z N )=C N , a l[\l-z l \ 2a H \z -z k \ 2 , Zl = e i6 '. 

1=1 l<j<k<N 
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The weight function has an algebraic singularity at z — e 10 = 1 corresponding in 
the log-gas picture of J3J to an impurity charge of a, so the ensemble is termed the 
spectrum singularity case. The probability that no eigenvalues e 1 ® have phases in 
the interval (—0, 9) is given by 



(4) E 2 (0;(-9,9y,\l-z\ 2a ;N) 

dz\ f dzN 



■p(zi, . . . , z N ). 



This problem can also be equivalently expressed in terms of a problem with its 
spectrum on the real line, where the weight function is now a Cauchy one, and the 
above gap probability is E 2 (0; (— oo, — s) U (s, oo); (1 + A 2 ) -JV ~ a ; N). 

We wish to focus here on the specific example of a double interval J = (— x,x) 
for the bulk scaling limit of the Spectrum Singularity N — ► oo. The gap 

probability is known to be expressed as a Fredholm determinant 

(5) £f s (0;(-ar,x);o)=det(l-Kj), 

where the integral operator IK has the kernel 

6 K(x, y) = Virxy/¥i] r , 

2(x - y) 

with a parameter a G C with Re(a) > —1/2 and the density of eigenvalues is 
p = 1. The Tracy and Widom theory for the Fredholm determinant forms of gap 
probabilities ^B] was employed in Ej , and this probability was evaluated as 

(7) Ei s (0;(-x,x);a) = exp ( [ ' "* dy CTl< '' " 



h y 

in terms of <J\{r) = —2xR(x, x) with r = 2x, which was shown to satisfy the 
ordinary differential equation 



(8) {ra'if +A[-a 2 -a 1 +ra' 1 ]{(a' 1 ) 2 - a-^fcP + ax 
subject to the boundary condition 

(9) o\ (r) 



= 0, 



2(2a + 3) (2a + 1) 16(2a + 5) (2a + 3) (2a + 1) 

r 4a+2 / a + 1 \ r 6a+3 



where 

(10) C a 



4^T(a + y 2 )T(a+r 2 )- 



When a — this equation reduces to case of the bulk scaling limit first found in 
|13| . and is then a special case of the Jimbo-Miwa-Okamoto a- function form of Py ■ 
Equation (jHJ was solved in terms of Painleve's fifth transcendent ^Hj for general a 
with the parameters 

(11) a = V 32 (l - 2af, - -Y 32 (l - 2a) 2 , 7 = 0, S = -2 . 
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However there is a puzzling aspect to this result, which also appears in the cases 
of double intervals with other weights such as the Gaussian, symmetric Jacobi [2U] 
or Cauchy ones. Expressed as a second-order second-degree ODE JSJ has rad- 
ical terms, namely the square-root term, or alternatively is a quartic second-order 
ordinary differential equation in polynomial form. Consequently, the ordinary dif- 
ferential equation (JHJ is not of the Jimbo-Miwa-Okamoto c-form of Py for a ^ 0, as 
is the case for the single intervals. On this point we remark there is now a large body 
of works [7| demonstrating that gap probabilities and their generalised aver- 
ages for unitary random matrix ensembles with classical weights can be evaluated 
in terms of a single t- function for one of the Painleve systems. Consequently the 
logarithmic derivatives of these averages satisfy a Jimbo-Miwa-Okamoto u-form for 
the appropriate system, which are generically only second-order second-degree or- 
dinary differential equations (see the classification of such ODEs with the Painleve 
property by Cosgrove and Scoufis 1 ) . Therefore it was not clear how the spectrum 
singularity result fitted into this broader scheme. We provide an answer to this 
question in this work - there is a natural explanation in the Okamoto theory for 
the Painleve transcendents Pm and Pup . 

Firstly we recognise that the special case (jl 1|> of Py is one that degenerates to 
Pin using the transformations of Gromak |11U12| - for example this can be achieved 
with the parameters (v\, 112) = (a — a — Yi) or am = 1 — 2a, /3m = 1 + 2a, 7m = 
1) £>m = — 1. Equation (jSJ arises quite naturally in the Painleve III system, as will 
be apparent from Okamoto's theory for Pm |15|. The Pm differential equation is 

(12) q" = - (q'f -a' + -{aq 2 + f3) + W + -, 

q t t q 

and can be generated from the Hamiltonian system 

(13) tH m = 2q 2 p 2 - [2 Voo tq 2 + (2vi + l)q - 2 Vo t] p + ^{Vi + v 2 )tq. 

One then recovers the standard Pm ODE for q(t) with the parameter identifications 

(14) a = -47700^2, (3 = 4?7 (wi + 1), 7 = 4r^, 5 = -Ar, 2 . 

Here rjoo , rj are arbitrary parameters which control the scaling of the independent 
and dependent variables, so they are usually fixed at some nominal value (unity). 
Now let us define the auxiliary Hamiltonian 

(15) h = tH ni + V 8 (2vx + l) 2 , 
and examine the time evolution of this. 

Theorem 1 (Proposition 1.9, |15p. The auxiliary Hamiltonian h(t) for Pm , as 
specified by l\15\) with parameters v\ , t>2 satisfies the ordinary differential equation 

(16) (th") 2 = [2(h - th')} |4(«') 2 + 16r/ ?7oo [2(h - th')} 

- 16?7o?7o e(w2 - Vi - l)y/2(h - th') - 16%77oo(w 2 - + V2)}, 

with an arbitrary sign e = ±1. 
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Proof: One can verify that the canonical variables are given by 

th" 



Ar] p = hi - e 

h'-, . 

y/&{h - th' 

2r\ x q — 



y/8(h - th') 
th" 



v 2 -Y 2 - ey/2(h-th') 

with e = ±1 provided that h(t) is not a singular solution of (|16|l . Employing these 
two relations in 

8(h-th') = (4qp-2v 1 - l) 2 , 
one arrives at the stated ordinary differential equation l|16|l . □ 

Remark 1. Okamoto gives a variant of this (modulo typographical mistakes) for 
e = 1 where it is rearranged and squared to render it polynomial. 

Remark 2. A more restricted version of this can be found using the Ablowitz-Fokas 
method of Riccati transformations |2J. Thus by extending their transformation 
(Theorem 4.3) relating the solution v{z) of a Pni equation with general parameters 
to that of a second-order quadratic ODE with solution </>(i), then defining a new 
variable w(t) by 

(jj, 2 = A(tw' -w)+r), 
one finds w(t) satisfies a particular form of (|16|l . 

By making scale changes to dependent and independent variables JSJ can be 
brought into correspondence to l|16|) and one solves for the parameters, 

(17) vi = -ea-y 2 , v 2 = ea + ]4 

(18) r] = 7]oo = 1, r = Ait 

(19) CTi (r) = 2h(t) - a 2 = 2tH m (t). 

Thus tTi (r) jr is proportional to the Hamiltonian for the Pm system, with parame- 
ters as specified. Furthermore, introducing the r-function Tni(i) for the Pm system 
by the requirement 

(20) H m (t) =: j t logmi(t), 
it follows from JTHJl, JSJl and Q that 

(21) B 2 ss (0;(-a;, a :);a)=r I 2 II (^) 

So far everything is in terms of the Pni system, but now we will find that 
conversion to the Pup system yields the desired r-function representation. The 
Pur system {s, 5m' , pup , ffur} is entirely equivalent to the Pni system but differs 
in a number of salient features, one of which is that the transcendent qiw(s) has 
movable double poles in contrast to the single poles for qm(t). 

The Hamiltonian for the Pm' system is given by |15| 

(22) sH uv (s) = qurPur - (ftn- + Uigup - s)puv + l A{vi + v 2 )quv, 
and the associated r-function is given by 

(23) H lIV =: -^-logT lIV {s). 

as 



From the Hamiltonian the cr-function ami is defined by 
(24) cr nv {s) 



(sHn V )(s/A) - - v 2 ) + js. 



The function anv(s) for general v = (vi,v 2 ) satisfies the cr-form second-order 
second-degree ordinary differential equation 

(25) (sa uv ) 2 - v lV2 (4w) 2 + °m>{^'m> - l)(^iH' - stria,) - " = 0. 

Lemma 1. The Hamiltonian for the Pm system is related to that of the Pm' system 
by 



= sHav(s) +sH uv (s) 

V V 

r a (v) = («i+i,«2-i), 



(26) tff m (t) 

w/iere 
(27) 
and 
(28) 

Proof. Using the mapping relating the two systems |15| 

(29) tq IU (t)=quv(s) 
1 



T 2 (v) 



t 2 = s. 



-Pm(t) =Piw(s) 



(30) 

(31) tH m (t) = 2sH U v{s) - quv(s) PlIV (s), 

we notice that the right-hand side of the last equality can be written 

(32) 2sH uv (s) - quvpuv = sH m ,(s) + {sH uv (s) - quvpuv) 

(33) =sH U v(s) + T 2 (sH m ,(s)) 



(34) 



up 



sHi 



(V1.V2) 



up 



('01 + l,l?2 — 1) 



where the Schlesinger transformation T 2 acts on the parameters according to l|27|) 
and upon the Hamiltonian according to T2(s-ffni') = s ^iip I v ^t 2 (v) ( see Table QJ. 
Thus we have ©. □ 

Proposition 1. The gap probability is the product of two t -functions for the 

system Pup 

(35) 



(-ea+y 2 ,ea-y 2 ) 



E$ s (0; (-x,x);a) = t iiv (-^(ttx) 2 ) m v (-yi(wx) 2 ) 

(-ea-Y 2 .ea+y 2 ) 

fore = ±1. The relevant solutions to i25\) appearing here are those with v = (/!, — /i) 
and /1 = — ea — y^ or ji = — ea + ^ satisfying the boundary condition 



(36) criii'(s) 



2 



+ " ( 1 

-a 2 



s 
4 



1 



2/i - 3 



2( M -2) 16(/i-3)(/x-2)( M -l) 
s) 2 ~ 2 ^ A 2/1-3 



/i — 1 



1 



2(M-2) 2 



where 
(37) 



C ^ i-^ri-n + 2)r(- M + 1) 



There is another known product form of (O U |H| which states 



(38) El s (0; {-x, x);a) = E™(0; (0, ttV); a - V 2 )E^(0; (0, irV); a + 



tiHE 



HE/ 



where the hard edge gap probability has the evaluation 
(39) 



E™(0;(0,X);a) = e- x ^mv(y 4 X) 



exp 



I d i am ' {s) ) 



v—(a,a) 



again in terms of r-functions for the Pup system. The particular solution of 1|25|) 
arising in the hard edge case satisfies the boundary condition 



(40) anv(s) ~ C a s a+1 

s^0+ 



1 - 



2a + 3 



2(a + 2) 16(a + 3)(a + 2)(a + 1)' 



„2a+2 



c. 



a a+l 



1 - 



2a + 3 
2(a + 2) 2 



„3a+3 



(1 + ...) + 



where 
(41) 



Ca 



2 2a + 2 r(a + 2)r(a+ 1)' 
Combining these two product forms we have the general identity 



(42) e- 2A r„p(X) 



( a -y 2 ,a-y 2 ) 



mv(x) 



(a+y 2 ,a+y 2 ) 



mv(-X) 



mv(-x) 



(-ea-y 2 ,ea+y 2 ) 

or an additive relation in terms of the Hamiltonian functions 



(~ea+y 2 ,ea-y 2 ) 



(43) -2s + sHiw(s) 



(a-y 2 ,a-y 2 ) 
= sHuv(-s) 



sH uv (s) 



(a+y 2 ,a+y 2 ) 

sHm>(—s) 



(-ea+y 2 ,ea-y 2 ) 



(-ea-y 2 ,ea+y 2 ) 

A direct derivation of the latter relation can be deduced from the actions of the 
reflection operators S2 and s_S2 which have the actions S2 '■ V2 ► —V2, s i— ► —s and 
s_S2 : vi — v\, s i — ► — s. Thus from Table 1 we see that 



(44) sH lw {-s) 

(45) sif nF (-s 



(— Vl,V 2 ) 



s-s 2 sHuv(s) 



(vi,v 2 ) 



(vt,v 2 ) 



{vi,v 2 ) 



sH uv (s) 



(vi,v 2 ) 



S, 



= S2sH uv (s) 

(vi, — v 2 ) 

and 143(1 follows by adding two instances of either of the above two relations. 

There is an interesting special case when the r-functions are classical solutions 
and this was found to occur in the studies ^3 d for a e Z>o + }/2- F° r the 
r-functions with parameters on the diagonal (a = n — y^, n £ Z>o) 



(46) 



thf(X) 



(a+y 2 , a +y 2 ) 



= dct 



I 3 - k (2y/X) 



j,k— 0,...,n — 1 



whereas on the cross-diagonal 



(47) mi'PO =e x Aet Jj_ k (2VX) 

(-a-y 2 ,a+y 2 ) L 



j,k— 0,...,n — 1 



where J v {z) 1 Iv(z) are the Bessel function and modified Bessel function respectively. 
Naturally the identity l|42(l is satisfied for these two sets of classical solution. 

In conclusion our results raise the natural question whether a similar evaluation 
in terms of r-function products as in (|35|l should exist for the gap probabilities 
defined on symmetric double intervals with the finite N ensembles of Gaussian, 
symmetric Jacobi and Cauchy weights. These are known to be solved in terms 
of Py and Pvi transcendents and there are analogues of the higher degree ODE 
©. Moreover the analogue of the product formula 138fl is known for each of these 
cases @], The next question is whether there is a corresponding Hamiltonian 
theory underlying these ODEs and how it relates to the known Hamiltonian theory. 
Isolated examples of higher degree second-order ODEs (quadratic and quartic in the 
second derivative) are known in the literature through many differing approaches 
|17j . 521 j E3 , EE] , and one may suspect that these too appear in a random matrix 
context. 
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